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The present study deals with an empirical investigation between CO, emissions, energy intensity, 
economic growth and globalization using annual data over the period of 1970-2010 for Turkish economy. 
We applied unit root test and cointegration approach in the presence of structural breaks. The direction 
of causality between the variables is investigated by applying the VECM Granger causality approach. Our 
results confirmed the existence of cointegration between the series. The empirical evidence reported that 
energy intensity and economic growth (globalization) increase (condense) COz emissions. The results 
also validated the presence of environmental Kuznets curve (EKC). The causality analysis shows 
bidirectional causality between economic growth and CO, emissions. This implies that economic growth 
can be boosted at the cost of environment. 
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1. Introduction 


Turkey has experienced a significant rise in economic growth, 
energy consumption and carbon emissions during the last two 
decades. Turkey is a candidate for full membership of European 
Union (EU) and therefore is likely to face significant pressures from 
EU during negotiations to introduce its national plan on climate 
change and global warming along with specific emissions targets [1]. 
Turkey is one of the important countries which have a high carbon 
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emissions and economic growth in the world. The reports of World 
Bank and UNDP indicate that CO emissions would rise more than 
six-fold by the end of 2025 rather than 1990s, so it is a great 
challenge for Turkey to achieve both the targets of high economic 
growth and less CO, emissions at the same time. 

The present study contributes in energy economics by four 
ways: (i), we augmented the CO2 emissions function by incorpor- 
ating globalization as a potential determinant of energy intensity, 
economic growth and CO2 emissions; (ii) Zivot-Andrews [2] unit 
root test has been applied in determining integrating order of the 
variables; (iii) Gregory-Hansen structural break cointegration test 
is used to examine the robustness of long run relationship 
between the variables; and (iv) direction of causal relation is 
investigated by applying the VECM Granger causality test. Our 
findings confirm the existence of long run relationship between 
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economic growth, energy intensity, globalization and CO2 emis- 
sions. We find that the EKC is validated in case of Turkey. More- 
over, energy intensity is a, major contributor to CO, emissions, and 
globalization improves the environmental quality. The feedback 
effect exists between economic growth and CO, emissions. Energy 
intensity and globalization Granger cause CO2 emissions. 

The rest of the paper is organized as follows: Section 2 presents 
review of literature; Section 3 provides data information, model- 
ing and estimation strategy; result interpretations are in Section 4 
and Section 5 deals with conclusion and policy implications. 


2. Literature review 


In 1991, Grossman and Krueger started the debate of environ- 
mental Kuznets curve (EKC) which explained the relationship 
between environmental pollution and economic growth i.e. 
inverted U-shaped relationship. Later on, a series of debate has 
started by investigating the relationship between environmental 
pollution and economic development. Johansson and Kristrém [3] 
noted that the literature on EKC is not enough and this topic needs 
more empirical investigation. But Stern [4] argued the issues of 
EKC should be revisited by using new models and new decom- 
positions with different panels and time series data sets. Similarly, 
Wagner [5] pointed out that the data on per capita CO emissions 
and per capita GDP are not stationary in time series framework 
and this problem is not sufficiently addressed in literature. There- 
fore, many dimensions of EKC are available for further empirical 
investigation. 

Existing literature provides two strands of relationship 
between energy consumption and energy emissions i.e. economic 
growth and energy consumption and, economic growth and CO2 
emissions in case of Turkey (see [6], for literature survey on 
energy-growth nexus). For example, Altinaya and Karagol [7] 
investigated the direction of causality between energy consump- 
tion and economic growth. They applied unit root test to examine 
stationarity properties of the variables. The Hsiao Granger caus- 
ality was applied using time series data over the period 1950- 
2000. Their empirical exercise reported the neutral effect between 
economic growth and energy consumption. Lise and Montfort [8] 
probed the relationship between gross domestic product and 
energy consumption using annual data for the period of 1970- 
2003. The Granger causality analysis found unidirectional causality 
running from economic growth to energy consumption. 

Soytas and Sari [9] analyzed the relationship between energy 
consumption in industrial and manufacturing sectors using multi- 
variate model by incorporating capital and labor in production 


Table 1 
Summary of studies on Turkey. 


function. Their results indicated cointegration between the vari- 
ables for long run relationship. The results of vector error correc- 
tion (VECM) model reveal that there is unidirectional causality 
running from energy consumption to manufacturing GDP. Further- 
more, the results of variance decomposition and generalized 
impulse response analysis confirmed that energy consumption is 
an important factor of manufacturing GDP. This implies that 
utilization of energy saving modes and energy efficiency technol- 
ogy may enhance manufacturing production in Turkey. Similarly, 
Jober and Karanfil [10] reinvestigated the relationship between 
energy consumption and economic growth at aggregate level and 
at sectoral industry level. Their results reported that there is no 
causality between both variables at aggregate level as well as 
sectoral level. 

Erdal et al. [11] used the data over the period of 1970-2006 to 
reexamine the relationship between energy consumption and real 
GNP. They applied augmented Dickey-Fuller (ADF) and Philips- 
Perron (PP) unit root tests to test stationarity properties of both 
variables and Johansen cointegration for long run as well as Granger 
causality test for pair-wise causality. The empirical exercise reported 
the cointegration between energy consumption and real GNP. The 
causality analysis revealed feedback effect implying that economic 
growth and energy consumption are interdependent. This suggests 
that any negative energy shock will put negative effect on economic 
growth of Turkey. Kaplan et al. [12] reexamined the causal relation- 
ship between economic growth and energy consumption over the 
period 1971-2006 using supply and demand side models. Their 
results found long run relationship as cointegration is found between 
the series. The causality analysis found feedback effect between 
economic growth and energy consumption. This shows that for 
achieving high level of economic growth more energy is needed 
and supply of energy further enhances economic growth, implying 
that any shock which occurred in supply of energy puts a negative 
impact on economic growth. 

The second strand deals with the relationship between eco- 
nomic growth and CO, emissions such as, Akbostanci et al. [13] 
tested the direction of causality between income and environ- 
mental degradation using various stages of economic development 
using PMj) and SO, measures of environmental degradation. They 
used data of 58 provinces of Turkey over the period 1968-2003. 
Their empirical results unveiled that CO emissions and income 
have long run relationship but inverted U-shaped relationship is 
found when they SO, and PMjo are used as measures of environ- 
mental degradation. The results do not support EKC hypothesis 
based on income and environmental degradation nexus. 

Halicioglu [14] augmented CO, emissions function by incorpor- 
ating trade to investigate the causal relationship between income, 


No. Authors Time period Variables Cointegration EKC hypothesis 

1 Altinaya and Karagol [7] 1950-2000 GDP per capita, energy consumption - - 

2 Lisea [19] 1980-2003 GNP growth, CO, emissions - - 

3: Lise and Montfort [8] 1970-2003 GDP per capita, energy consumption Yes - 

4. Soytas and Sari [9] 1968-2002 Energy consumption, industrial value-added Yes - 

5 Jober and Karanfil [10] 1960-2003 GNP growth, energy consumption - - 

6 Erdal et al. [11] 1970-2006 Energy consumption, GNP per capita Yes - 

7. Akbostanci et al. [13] 1992-2001, GDP per capita, CO2 emissions Yes EKC does not exist 
1968-2003 

8. Halicioglu [14] 1960-2005 Income, CO, emissions, energy consumption, trade openness Yes - 

9. Soytas and Sari [15] 1960-2000 GNP, CO2 emission, energy consumption Yes - 

10. Ozturk and Acaravci [1] 1968-2005 GNP per capita, CO2 emissions, energy consumption, employment Yes EKC does not exist 

11. Jobert et al. [17] 1970-2008 GNP per capita, CO,emissions, energy consumption - EKC exists 

12. Kaplan et al. [12] 1971-2006 GNP per capita, energy consumption Yes - 

13. Joberta and Karanfil [18] 1971-2008 Energy consumption, CO2 emissions, real GDP per capita -= EKC exists 

14. Ozturk et al. [20] 1960-2006 GNP growth, energy consumption Yes - 
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CO2 emissions and energy consumption for the period 1960-2005. 
Halicioglu found cointegration by applying the ARDL bounds 
testing approach to cointegration. The results showed that GDP 
is highly significant among other variables of the model in 
explaining CO, emissions in case of Turkey. Soytas and Sari [15] 
reexamined the relationship between economic growth, CO2 
emissions and energy by incorporating capital formation and labor 
as potential determinants of economic growth and CO, emissions. 
Their results exposed that CO2 emissions Granger cause energy 
consumption but same and vice versa which implies that by 
reducing CO» emissions, Turkey may not forgo economic growth. 
Kaygusuz [16] investigated the electricity and energy demand 
functions and their empirical exercise found that rapid energy 
consumption and energy production are linked with environmen- 
tal issues at national level as a rise in energy consumption 
(electricity consumption) increases CO2 emissions. 

Ozturk and Acaravci [1] reinvestigated the cointegration and 
causality between economic growth, CO emissions and energy 
consumption by incorporating employment using time series data 
over the period 1968-2005. Their results indicated the existence of 
cointegration between the variables and found that income 
elasticity of CO, emissions is inelastic (-0.606) but income 
elasticity of energy consumptions is more elastic (1.375). Their 
analysis could not provide the empirical validation of the EKC 
hypothesis. The causality analysis found neutral effect between 
energy consumption and economic growth, economic growth and 
CO- emissions and, energy consumption and CO, emissions. This 
implies that the adoption of energy conservation has no adverse 
effect on growth rate of real GDP. Jobert et al. [17] probed the 
relationship between economic growth, CO2 emissions and energy 
consumption by applying Bayesian empirical model. The study 
used time series data of 50 countries including Turkey over the 
period of 1970-2008. Their empirical analysis reported that 
existence of the EKC is sensitive with respect to countries but 
EKC exists in case of Turkey. Joberta and Karanfil [18] used cross- 
country data including Turkey to test the validation of EKC and 
found the existence of EKC before 1980. The threshold level of 
income is rising which was reported as 10,000 in early 1980 and 
20,000 in 2008 (Table 1). 

The 21st century has increased the internationalization among 
the world economies and countries are more closely linked with 
each other economically, politically and culturally. Globalization 
which is based on economics facilitates and helps in the promo- 
tion of division of labor and increases the comparative advantage 
of different nations. Globalization improves the total factor pro- 
ductivity by increasing trade activity but also boosts economic 
activity via foreign direct investment and transfer of advanced 
technology from developed countries to developing nations. 
Globalization also provides investment opportunities including 
foreign direct investment and develops the financial markets. 
Globalization directly enhances trade and then economic growth 
while indirectly, it promotes investment opportunities not only in 
the form of domestic investment but also in foreign investment, 
which not only influences energy demand but also influences the 
environment. 

Various researchers have used different measures of globaliza- 
tion to examine its impact on environmental degradation. For 
instance, Grossman and Krueger [21] investigated the environ- 
mental impact of Northern America Free Trade Agreement 
(NAFTA) on environment. They reported that trade openness 
(globalization) affects environmental degradation via scale effect 
keeping composition effect and technique effect constant. Simi- 
larly, Dinda [22] claimed that environmental degradation increases 
as scale effect dominates the composition effect and technique 
effect and same conclusion is drawn by Shahbaz et al. [23] that 
trade openness declines CO2 emissions in case of Pakistan. On the 


contrary, Wheeler [24] noted that globalization reduces environ- 
mental degradation due to the investment in energy-efficient 
technologies for production. Copeland and Taylor [25] reported 
that globalization facilitates transfer pollution intensive technol- 
ogy to countries where environmental regulation are weak (in 
developing economies). In such circumstances, developed coun- 
tries attain benefits from trade openness at the cost of environ- 
ment in developing economies. Copeland and Taylor [26] pointed 
that trade depends upon the relative abundance of factor endow- 
ment in each country and therefore, comparative advantage of 
trade also affects environmental quality depending upon trade and 
environmental policy in the country. Birdsall and Wheeler [27], 
Lee and Roland-Host, [28], Jones and Rodolfo [29] opined that 
environmental degradation is not the main cause of trade open- 
ness. Similarly, Antweiler et al. [30] and Liddle [31] pointed out 
that trade openness improves environmental quality via technique 
effect. Environmental regulations become strict as income 
increases and the adoption of energy-efficient technologies are 
encouraged to save environment from degradation. In case of 
China, Dean [32] reported that trade openness deteriorates envir- 
onmental quality via improved terms of trade, however, rise in 
income saves environment from degradation. Magani [33] used 
data of 63 developed and developing economies to examine the 
effect of trade openness on energy emissions. The results showed 
that a 0.58% carbon emission is linked with a 1% increase in trade. 
Similarly, McAusland [34] reported that trade affects environment 
significantly and same view is confirmed by Frankel [35]. 


3. The data, modeling and estimation strategy 
3.1. The data and modeling 


We have used data of energy intensity per capita, CO emis- 
sions per capita, real GDP per capita and globalization index to 
probe the existence of environmental Kuznets curve (EKC) in case 
of Turkey. The data on energy consumption (kilotons of oil 
equivalent), CO2 emissions (metric tons) and real GDP (Turkish 
currency) has been attained from world development indicators 
(CD-OM, 2012). The series population is used to convert all series 
into per capita. The data on KOF globalization index is borrowed 
from Dreher [36]. The globalization index is constructed from 
three sub-indices (social, economic and political globalization).? 
The study covers the period of 1970-2010. The general functional 
form of our model is given in the as following equation: 


ACr =f(Et, Ye, Y?, Ge) (1) 


We have transformed all the variables into natural logarithm 
following Shahbaz et al. [23,37]. The empirical form of our model 
is constructed as follows: 


InC; = ay + agInE; + aylny; 4 ayaInY? + aglnGe + py (2) 


where InC; is natural log of CO2 emissions per capita, natural log 
energy intensity per capita is indicated by InE;, InY; (InY?) is the 
natural log of real GDP per capita (square of real GDP per capita) 
and InG; is for natural log of KOF index of globalization. y, is the 
error term assumed to be having normal distribution with zero 
mean and predictable variance. We expect that impact of energy 
consumption on CO, emissions and ag>0. The relationship 
between economic growth and CO2 emissions inverted U-shaped 
if ay >0 and ayz <0 otherwise U-shaped if ay <0 and ay > 0. 
ac < 0 if energy-efficient technology via foreign direct investment 
and trade is encouraged for domestic production otherwise ag > 0. 


? See in details http://globalization.kof.ethz.ch/. 
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3.2. Zivot-Andrews unit root test 


Numerous unit root tests are available to test the stationarity 
properties of the variables including ADF by Dickey and Fuller [38], 
P-P by Philips and Perron [39], KPSS by Kwiatkowski et al. [40], 
DF-GLS by Elliott et al. [41] and Ng-—Perron by Ng-Perron [42]. 
These tests provide biased and spurious results due to non- 
availability of information about structural break points in series. 
In doing so, Zivot-Andrews [2] developed three models to test the 
stationarity properties of the variables in the presence of struc- 
tural break point in the series: (i) this model allows a one-time 
change in variables at level form, (ii) this model permits a one- 
time change in the slope of the trend component i.e. function and 
(iii) model has one-time change both in intercept and trend 
function of the variables to be used for empirical analysis. Zivot- 
Andrews [2] followed three models to validate the hypothesis of 
one-time structural break in the series as follows: 


k 
AX; =A + AX,_1 + bt + cDU; + > djAxt-j + He (2) 
j=1 
k 
AX, = b+ Dxp_4 +ct+ bDT; + bs dj AX; + He (3) 
j=l 
k 
AXt = C+ CXt-1 + Ct + dDU; + dDT; + p2 djAX;-j + Ht (4) 


j=l 


where dummy variable is indicated by DU; showing mean shift 
occurred at each point with time break while trend shift variables 
is show by DT;.? So, 


1 if t>TB t-TB if t>TB 
DU; = 


0 if t<TB 


o if tetas 4 Dur={ 


The null hypothesis of unit root break date is c=0 which 
indicates that series is not stationary with a drift not having 
information about structural break point while c< 0 hypothesis 
implies that the variable is found to be trend-stationary with one 
unknown time break. Zivot-Andrews unit root test fixes all points 
as potential for possible time break and estimates through 
regression for all possible break points successively. Then, this 
unit root test selects that time break which decreases one-sided 
t-statistic to test ¢(=c-1)=1. Zivot-Andrews report that in the 
presence of end points, asymptotic distribution of the statistics is 
diverged to infinity point. It is necessary to choose a region where 
end points of sample period are excluded. Further, Zivot-Andrews 
suggested the trimming regions i.e. (0.15T, 0.85T). 


3.3. The ARDL cointegration 


We employ the autoregressive distributed lag (ARDL) bounds 
testing approach to cointegration developed by Pesaran et al. [44] 
to explore the existence of long run relationship between eco- 
nomic growth, energy intensity, globalization and CO, emissions 
in the presence of structural break. This approach has multiple 
econometric advantages. The bounds testing approach is applic- 
able irrespective of whether variables are (0) or I(1). Moreover, a 
dynamic unrestricted error correction model (UECM) can be 
derived from the ARDL bounds testing through a simple linear 
transformation. The UECM integrates the short run dynamics with 
the long run equilibrium without losing any long run information. 
The UECM is expressed as follows: 


AlnC; = ay +arT + aclnC;_4 + aglnE,_1 + aylnY;_1 


3 We used model-4 for empirical estimations following Sen [43]. 


p 
+ay2InY? , + aclnG;_1 + D> a, AlnCy_; 
i=1 
p r s 2 
+ > aj AInE;_; + > ay AlnY s-k + >. ajAlnY$ | 
0 


j=l k= 1=0 


J 
+ >, amAINGt_m + apD, + Hat (9) 


m=0 
AlInE; = ay +arT + aclnC;_4 + aglnE;_1 + aylnY;_1 


p 
+ay2 InY? , + aclnG,_1 + >, aj AlnE;_; 
tat 


q r S 
+ ¥ ajAlnCj+ X akAlnY k+ X aAlnY?, 
j=1 k=0 l=0 


t 
+ È amAlNGi-m + apD2 + po (10) 


m=0 
Alny; = ay +arT + aclnCy_1 + aglnE;_4 + aylnY 1 


Pp 
tay lnY?_, +aclnGei + È ANY, ; 
i=l 


q r $ 
+ È pAlnC,;+ Y pAlnEe+ X AAlnY? 
j=1 k=0 1=0 
t 
+ È BmANGm + BpD3 + Hzr (11) 
0 


m= 


AlnY? =a, +arT + aclnC;_4 + aglnE,_1 + aylnY;_1 


p 
+ay21nY? ; + aclnGy_1 + > 9;AlnY?_, 
i=1 


q r Ss 
+ > djAlnC,_; + > 9, AINE;_~ + > 9, AlnY;_) 
k=0 1=0 


j=0 


t 
+ È dmAING_1 + 9pD4 + pat (12) 


m=0 
AlnG; = ay +arT + aclnC;_4 + aglnE;_4 + aylnY 1 


Pp 
taylnY?_, + aclnNG + È p)AlnG,; 
i=1 


q r 
+ > pjAlnG,_; + > pyAlnE;_, 
j=0 k=0 


S t 
+ 2 pıAlnY -1 + z PmAINY? n + OpDs + fist (13) 
=0 m= 
where A is the first difference operator, D is dummy for structural 
break point and „is error term assumed to be independently and 
identically distributed. The optimal lag structure of the first 
differenced regression is selected by the Akaike information 
criteria (AIC). Pesaran et al. [44] suggests F-test for joint signifi- 
cance of the coefficients of the lagged level of variables. For 
example, the null hypothesis of no long run relationship between 
the variables is Ho : ac = ag = ay = ayz = ag =0 against the alter- 
native hypothesis of cointegration Hg : actagtay+#ay2#ac¢+0. 
Pesaran et al. [44] computed two sets of critical value (lower and 
upper critical bounds) for a given significance level. Lower critical 
bound is applied if the regressors are (0) and the upper critical 
bound is used for I(1). If the F-statistic exceeds the upper critical 
value, we conclude in favor of a long run relationship. If the F- 
statistic falls below the lower critical bound, we cannot reject the 
null hypothesis of no cointegration. However, if the F-statistic lies 
between the lower and upper critical bounds, inference would be 
inconclusive. When the order of integration of all the series is 
known to be I(1) then decision is made based on the upper critical 
bound. Similarly, if all the series are [(0), then the decision is made 
based on the lower critical bound. To check 
the robustness of the ARDL model, we apply diagnostic tests. 
The diagnostics tests check normality of error term, serial 
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Table 2 
Descriptive statistics and correlation matrix. 


Variables Inc; InY; InE; InG; 
Mean 0.8961 6.7912 6.7916 3.8872 
Median 0.9773 6.8019 6.8225 3.8954 
Maximum 1.5838 7.2413 7.2224 4.2481 
Minimum 0.2495 6.3641 6.2651 3.5562 
Std. dev. 0.3441 0.2609 0.2560 0.2368 
Skewness -0.0055 0.1553 -0.1652 -0.1119 
Kurtosis 2.0783 1.8726 1.8994 1.3856 
Jarque-Bera 1.4160 2.2789 2.2007 4.4271 
Probability 0.4926 0.3199 0.3327 0.1093 
Inc; 1.0000 

InY; 0.6083 1.0000 

InE; 0.7010 0.6901 1.0000 

InG; -0.3353 -0.0666 -0.0710 1.0000 


correlation, autoregressive conditional heteroscedasticity, white 
heteroscedasticity and the functional form of empirical model. 


3.4. The VECM Granger causality 


After examining the long run relationship between the vari- 
ables, we use the Granger causality test to determine the causality 
between the variables. In case of cointegration between the series, 
the vector error correction method (VECM) can be developed as 
followings: 


AlnC; b B11,1B121B13,1B14,1B12,1 AlnC;- 
AInE; Boi 1B 22,1 B23,1B241B 251 AInE;_1 
AlnY; | = “i + | B31,1B32,1B33,1B34,1B35,1 | x | AlnY;4 
AlnY? i Bay 1B421B43,1B44,1B45,1 AlnY? , 
AlnG; : B51,1B52,1 B531 B54,1B55,1 AlnG;_1 
B11,mB12,mB13,mB14,mB15,m 
B21 ,mB22,mB23,mB24,mB25,m 
+... + B31 mB32,mB33,mB34,mB35,m 
Bai mB42,mB43,mB44,mB45,m 
B51 ,mB52,mB53,mB54,mB55,m 
AlnCy_1 a Hat 
AInE;_4 62 Hat 
x | AlmYi-1 | + | 63 | x (ECM,-1)+ | Fat (14) 
AInY;., C4 Hat 
AlnG;_1 os Hst 


where difference operator is (1-L) and ECM;-, is the lagged error 
correction term, generated from the long run association. The long 
run causality is found by the significance of coefficient of lagged 
error correction term using t-test statistic. The existence of a 
significant relationship in first differences of the variables provides 
evidence on the direction of short run causality. The joint 
x Statistic for the first differenced lagged independent variables 
is used to test the direction of joint long-and-short runs causality 
between the variables. For example, B,2;+0V; shows that energy 
intensity Granger causes CO, emissions and energy intensity is 
Granger cause of CO, emissions if Bz, ;#0Vj;. 


4. Results discussion 


Descriptive statistics and correlation matrix are presented in 
Table 2. Based on Jarque—Bera test statistics, our results indicate 
that all the series are normally distributed having zero mean while 
variance is constant. This leads us to peruse for further analysis. 
The correlation matrix reveals a positive association between the 
underlying variables. For instance, economic growth is positively 


Table 3 
Zivot-Andrews structural break trended unit root test. 


Variable At level At 1st difference 
t-Statistic Time break t-Statistic Time break 
InC; -3.768 (1) 2001 -6.544 (1)* 1998 
InY; -4.210 (2) 1986 -6.559 (0)* 2004 
Iny? -4.108 (2) 1986 -6.684 (0)* 2004 
InE; -3.676 (1) 1986 -6.483 (0)* 2006 
InG; -2.927(0) 2005 -8.112 (1)* 1992 


* Significant at 1% level. Lag order is shown in parenthesis. 


correlated with CO, emissions and same is true between energy 
intensity and CO, emissions. A positive and high correlation is 
found between energy intensity and economic growth. Globaliza- 
tion is inversely correlated with CO, emissions, economic growth 
and energy intensity. 

There is a need to test the order of integration of the variables 
before applying the ARDL bounds testing to investigate long run 
relationship. Although, the ARDL bounds testing approach 
assumes that the variables must be stationary at I(0) or I(1) or 
1(0)/I(1). The process to compute F-statistic becomes invalid if any 
series is found to be stationary at I(2). So just to ensure that none 
of the variables is integrated beyond mentioned order of integra- 
tion, we have applied ADF unit root test. Our results of ADF and PP 
tests confined that none of the variables is stationary at level with 
intercept and trend. All the series are found to be integrated 
at I(1).4 The results of these may be biased and unreliable because 
both unit root tests ignore the role of structural break stemming in 
the series. The appropriate information about structural break 
arising in the series would be helpful to policy makers in 
articulating a comprehensive energy, economic, trade and envir- 
onmental policy to sustain long run economic growth. This issue 
has been resolved by applying Zivot-Andrews unit root test 
accommodating the information about single unknown structural 
break. 

The results shown in Table 3 reported that CO emissions, 
economic growth, energy intensity and globalization have unit 
root problem at level but all the series are stationary at 1st 
difference with intercept and trend. The same level of integrating 
properties of the variables intends to apply the ARDL bounds 
testing approach to cointegration to examine the long run rela- 
tionship between economic growth, energy intensity, globaliza- 
tion and CO, emissions over the period of 1971-2010 in case of 
Turkey. The two step ARDL procedure requires appropriate lag 
order of the variable to calculate F-statistic. We have used 
numerous lag length criterions and results are reported in 
Table 4. We lag selection is based on Akaike information criteria 
(AIC). Ltitkepohl [45] argued that AIC has superior power proper- 
ties for small sample data compared to other lag length criterions. 
The results reported in Table 4 noted that lag 2 is sufficient for 
such small sample data having only 40 observations (see third row 
of Table 4). 

The second step deals with the calculation of F-statistic to 
confirm whether cointegration between the variables i.e. eco- 
nomic growth, energy intensity, globalization and COzemissions 
exists or not. Table 4 presents the results of the ARDL bounds 
testing analysis. Our empirical evidence reveals that upper critical 
bound is less than our calculated F-statistic when we used CO; 
emissions, real GDP per capita (square of real GDP per capita) and 
energy intensity as dependent variables. Our F-statistics 11.656, 
8.941 (8.005) and 12.382 are statistically significant at 1(5)% level 


4 Results are available upon request from authors. 


M. Shahbaz et al. / Renewable and Sustainable Energy Reviews 25 (2013) 494-502 499 


Table 4 
Results of ARDL cointegration test. 


Estimated models Cc=f(¥t, Y?, Er, Ge) 


F-statistics 11.656* 8.941** 
Lag order 2.232622 2, 2,.2, 2,1 
Critical values 1% level 5% level 
Lower bounds 10.150 7135 
Upper bounds 11.130 7.980 
Diagnostic tests 

R? 0.9017 0.9997 
Adj-R? 0.7542 0.9994 
Durbin-Watson 2.0578 1.9480 


Ye =f(Ce, Y?, Et, Ge) 


Y? =f(Ct, Ye, Et, Ge) Er=f(Ce, Ye, Y?, Go) Gr=f(Ct, Ye, Y?, Ge) 


8.005** 12.382* 1.3342 

2, 2, 2, 2)-1 2,2; 1,21 Zi 2y Did) 2 
10% level* 

5.950 

6.680 

0.9997 0.9819 0.8018 
0.9994 0.9590 0.6037 
1.9572 1.8956 2.3233 


* Critical values bounds are from [47] with unrestricted intercept and unrestricted trend. 


* Significance at 1% level. 
** Significance at 5% level. 


Table 5 
Results of Johansen cointegration test. 


Hypothesis Trace statistic Maximum eigen value 
R=0 95.3434* 69.8188* 

Rs1 49.7529** 47.8561** 

Rs2 19.3135 29.7970 

Rs3 6.9000 15.4947 

Rs4 1.3126 3.8414 


* Significant at 1% level. 
** Significant at 5% level. 


of significance. This implies that we have four cointegrating 
vectors confirming long run relationship between economic 
growth, energy intensity, globalization and CO2 emissions over 
the period of 1971-2010. 

The robustness of long run results is investigated by applying 
Johansen and Juselies [46] cointegration approach and results are 
reported in Table 5. The results indicate two cointegrating vectors 
again confirming cointegration between the variables. 

The problem with results of the ARDL bounds testing devel- 
oped by Pesaran et al. [44] and Johansen and Juselies [46] is that 
they do not have information about structural break stemming in 
the series. Therefore in order to overcome this problem we have 
applied Gregory—Hansen [48] cointegration approach accommo- 
dating single structural break pointed out by Z-A unit root test. 
This test provides consistent and reliable empirical evidence as 
compared to other traditional cointegration tests [49]. Table 6 
reports the results of Gregory—Hansen cointegration and we find 
cointegrating single vector once we use economic growth, globa- 
lization and CO, emissions as forcing variables. This implies that 
cointegration is found in energy intensity equation after allowing 
for structural break in 1986. This break point is due to the usage of 
more coal instead of oil due to the oil crisis of 1970s. Overall, we 
find that long run results are robust. 

Table 7 deals with the long-run marginal impacts of economic 
growth, energy intensity and globalization on CO, emissions. The 
results expose that linear term of real income per capita has 
positive impact on CO emissions and whereas negative effect of 
square term of real income per capita on CO, emissions is reported 
which is statistically significant at 5% level of significance. This 
implies that inverted U-shaped relationship exists between real 
income per capita (square of real GDP per capita) and CO, 
emissions. The estimates of linear and nonlinear terms are 
7.3502 and -0.4336. This empirical exercise validates the existence 
of environmental Kuznets curve (EKC). This shows that a 1% 
increase in real income per capita is linked with 7.3502% increase 
in COz emissions and inverse effect of squared term of real income 


per capita indicates the delinking point of CO, emissions 
i.e.-0.4332, once an economy achieves threshold level of real 
income per capita. This justifies for the support of EKC which 
reveals that economic growth increases CO, emissions initially and 
improves the environmental quality once economy is achieves 
threshold level of income per capita. A positive relationship is 
found between energy consumption and CO2 emissions and it is 
statistically significant at 5% level. All else is remaining the same, a 
1% increase in energy consumption raises CO, emissions by 
0.7155% which shows that energy intensity is a major contributor 
to CO, emissions. Globalization has inverse impact on CO, emis- 
sions and is statistically significant at 10% level of significance. The 
results report that a 0.1950% decline in CO, emissions is due to 1% 
increase in globalization by keeping other things constant. 

The short run results are illustrated in lower part of Table 7. 
The results intend that linear and nonlinear terms of real GDP per 
capita have positive and negative signs (inverted-U shaped rela- 
tion) on CO, emissions and are statistically significant at 1% level 
of significance. The impact of energy intensity is positive on CO2 
emissions and it is statistically significant at 1% significance level. 
Globalization has inverse impact on CO, emissions at 10% level of 
significance. The coefficient of ECM;-1 has negative sign and 
significant at 10% level of significance. The significance of lagged 
error term corroborates the established long run association 
between the variables. Furthermore, the negative and significant 
value of ECM,- implies that any change in CO, emissions from 
short run towards long span of time is corrected by 27.54% every 
year. Sensitivity analysis indicates that short run model passes all 
diagnostic tests i.e. LM test for serial correlation, ARCH test, 
normality test of residual term, white heteroscedasticity and 
model specification successfully. The results are shown in lower 
segment of Table 7. It is found that short run model does not show 
any evidence of non-normality of residual term and implies that 
error term is normally distributed with zero mean and covariance. 
The serial correlation does not exist between error term and CO, 
emissions. There is no autoregressive conditional heteroscedasti- 
city and the same inference is drawn about white heteroscedas- 
ticity. The model is well specified proved by Ramsey RESET test. 
The stability of the ARDL bounds testing approach estimates is 
investigated by applying the CUSUM and CUSUMszq tests. The 
results are shown in Figs. 1 and 2. The plots of the CUSUM statistics 
are well within the critical bounds. 

The plots of the CUSUMsq test are not within the critical 
bounds. Furthermore, we apply Chow forecast test to examine 
the significance structural breaks in an economy for the period of 
2000-2010. In this study, F-statistic computed in Table 8 suggests 
that no significant structural break exists in case of Turkey during 
the sample period. The Chow forecast test is more reliable and 
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Table 6 
Gregory—-Hansen structural break cointegration test. 


Model Tc(C/Ye, Y?, Et, Gr) Ty(Yt/Cr, Y?, Et, Gr) 
ADF-test -3.0963 =3,8733 

Break year 2001 1986 

Prob. values 0.0028 0.0003 


Tye (Y? /Cr, Yt, Et, Gr) 


Te(Et/Ct, Yt, Y?, Gr) To(Gt/Ct, Yt, Y?, Et) 


-3.7339 -4,9861** -3.5329 
1986 1986 2005 
0.0004 0.0006 0.0008 


The ADF statistics show the Gregory-Hansen tests of cointegration with an endogenous break in the intercept. Critical values for the ADF test at 1%, 5% and 10% are -5.13, 


-4.61 and -4.34 respectively. 
** Significance at 5% level. 


Table 7 
Long-and-short run analysis. 


Dependent variable =In C, 


Long run results 


Variable Coefficient Std. error t-Statistic 
Constant -28.2399* 9.4667 -2.9830 
InY; 7.3502** 2.8005 2.6245 
Iny? -0.4336"* 0.2030 -2.1362 
InE; 0.7155" 0.3162 2.2627 
InG; =0,1950"" 0.1149 -1.6968 
Short run results 

Constant 0.0130 0.0104 1.2448 
InY; 10.3856* 3.3239 3.1245 
Iny? -0.7473* 0.2426 -3.0801 
InE; 0.7696* 0.1913 4.0223 
InG; -0.2952*** 0.1686 -1.7504 
ECM; -0.2754""* 0.1578 -1.7452 
Diagnostic tests 

Test F-statistic Prob. value 

x? NORMAL 3.6450 0.1616 

y SERIAL 0.9255 0.3434 

7° ARCH 0.9063 0.3476 

x? WHITE 0.7055 0.7110 

y RAMSEY 2.6867 0.1112 


* Significant at 1% level. 
** Significant at 5% level. 
**** Significant at 10% level. 


1980 1985 1990 1995 2000 2005 2010 


Fig. 1. Plot of cumulative sum of recursive residuals. The straight lines represent 
critical bounds at 5% significance level. 


preferable than graphs [50]. This confirms that the ARDL estimates 
are reliable and efficient. 

The presence of cointegration among the variables implies that 
causality relation must exist at least from one side. The directional 
relationship between energy intensity, economic growth, globali- 
zation and CO emissions will provide help in articulating com- 
prehensive policy to sustain economic growth by controlling 
environment from degradation and utilizing energy efficient 
technologies imported from advanced countries. We applied 
Granger causality test within the VECM framework to detect the 
causality between the variables. Table 9 reports the results of the 
VECM Granger causality analysis. The long run causality is cap- 
tured by a significant t-test on a negative coefficient of the lagged 
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-0.4 r 
1980 1985 1990 1995 2000 2005 2010 


Fig. 2. Plot of cumulative sum of squares of recursive residuals. The straight lines 
represent critical bounds at 5% significance level. 


Table 8 
Chow forecast test. 


Chow forecast test: forecast from 2000 to 2010 


0.2090 
0.0556 


1.5110 
13.7608 


F-statistic 
Log likelihood ratio 


Probability 
Probability 


error-correction term ECM;-1. The jointly significant LR test on the 
lagged explanatory variables shows short-run causality. 

Table 9 reveals that the estimates of ECM;-; are having negative 
signs and statistically significant in all the VECMs. The significance 
of lagged error term shows speed of adjustment from short run 
toward long run equilibrium path in the equation of energy 
consumption (-0.7311), economic growth (-0.4990, -0.4955) as 
well as CO2 emissions (-0.3801). This implies that the VECM 
equation of energy intensity has high speed of adjustment 
(-0.7311) as compared to economic growth (-0.4990, -0.4955) 
and CO2 emissions (—0.3801) the VECMs. 

The long run causality result reported that the feedback effect 
is found between economic growth and CO emissions. This 
implies that Turkey is achieving economic growth at the cost of 
environment. The bidirectional causality is found between energy 
intensity and CO emissions. This suggests adopting energy- 
efficient technology to enhance production which emits less CO2 
emissions. The feedback effect also exists between economic 
growth and energy intensity. 

This reveals that energy is an important stimulant like other 
factor of production and reduction in energy consumption would 
retard economic growth. This finding supports for energy explora- 
tion policies to sustain economic growth for long run. Globaliza- 
tion Granger causes economic growth, energy intensity and CO, 
emissions validating the globalization-led growth, globalization- 
led energy and globalization-led CO emissions hypotheses. 

The results of short run causality are very interesting. The 
neutral effect is found between CO emissions and globalization 
and, the same is true for energy intensity and globalization. 
Economic growth Granger causes CO2 emissions. The feedback 
effect exists between energy intensity and CO, emissions. There 
is no causality between energy intensity and economic growth. 
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Table 9 
VECM granger causality analysis. 


Dependent Direction of causality 
variable 
Short run Long run 
AlnC;-; AlnY;-1 AlnY2, AInE;-;  AlnG;-1 ECT; 
AlnC; - 3.0346*** 2.9838*** 14.2094* 0.3061 -0.3801""* — 
[0.0672] [0.0688] [0.0001] [0.7390] [-1.7435] 
AlnY; 0.7304 - 78.5723* 1.0836 2.8899*** -0.4990* 
[0.4913] - [0.0000] [0.3532] [0.0736] [-4.1581] 
AlnY? 0.6585 82.1544* - 0.9978  3.1674*** -0.4955* 
[0.5282] [0.0000] - [0.3824] [0.0587] [-4.2409] 
AlnE; 8.7686* 0.2828 0.2198 - 0.5230 -0.7311* 
[0.0012] [0.7558] [0.8041] - [0.5986] [-1.8901] 
AlnG; 1.2117 1.3052 1.1836 0.0566 - - 
[0.3139] [0.2883] [0.3221] [0.9450] 


* Significance at 1% level. 
** Significance at 5% level. 
** Significance at 10% level. 


The joint causality results are also reported in Table 9 validating 
our long-and-short runs findings. 


5. Conclusion and policy implications 


This paper probes the relationship between CO, emissions and 
economic growth by incorporating energy intensity and globaliza- 
tion as potential determinants of economic growth and CO, 
emissions in case of Turkey using annual data over the period of 
1970-2010. We have applied cointegration approaches to test the 
robustness of long run relationship between the variables is in the 
presence of structural breaks. The VECM Granger causality 
approach has been applied to examine the causal relationship 
between economic growth, energy intensity, globalization and CO3 
emissions. 

Our empirical exercise confirms the existence of cointegration 
in the presence of structural breaks in the series. Moreover, 
inverted U-shaped relationship is found between economic 
growth and CO emissions i.e. Environmental Kuznets curve 
(EKC). Thus, beyond a threshold level of real GDP per capita, any 
increase in real GDP per capita is likely to reduce the carbon 
emissions per capita in Turkey. 

Energy intensity increases CO, emissions and is a major 
contributor to energy emissions. Globalization seems to lower 
CO, emissions. We find feedback effect between economic growth 
and CO, emissions. The bidirectional causality is found between 
energy intensity and CO emissions and same is true for energy 
intensity and economic growth. Economic growth, energy inten- 
sity and CO2 emissions Granger cause globalization. 

Turkey is a candidate for full membership to the European 
Union (EU) and signed Kyoto Protocol to introduce its national 
plan on climate change and global warming along with specific 
emission targets and the associated abatement policies. Thus, 
numerous measures had been taken in last few years. However, 
these measures are not adequate for reducing environmental 
pollution without any sacrifices on the Turkish economic growth. 
To decrease carbon emissions and energy import related current 
account deficit of Turkey, the usage of alternative energy sources 
(renewable energy sources) like solar, wind, geothermal sources 
and bio-diesel fuel should be increased and green investment 
technologies should be supported. 

For future research, renewable and non-renewable energy 
sources of energy can be incorporated in neo-classical production 


501 


Joint long-and-short runs causality 


AlnC;-1, ECT;-1 AlnY;,-1, ECT,-1 AlnY2.,, ECT; AlnE,-1, ECT;-1 AlnG,-1, ECT;-1 


2.9935 2.9572** 14.7276" 2.4427" 

[0.0449] [0.0517] [0.0000] [0.0877] 
6.9671* = 58.0545* 6.0488* 6.2609* 
[0.0014] = [0.0000] [0.0029] [0.0024] 
7.0971* 58.9627* - 6.2364* 6.4747* 
[0.0012] [0.0000] - [0.0025] [0.0020] 
7.2650" 3.2025 3.2027" 2 4.4582** 
[0.0010] [0.0213] [0.0212] 2 [0.0114] 


to examine the relationship between energy consumption and 
economic growth following Shahbaz et al. [51] and Leitão [52] by 
incorporating the globalization. Globalization is a potential deter- 
minant of economic growth and energy consumption. There is 
need of empirical investigation of sectoral environmental Kuz- 
nets's curve in Turkey to improve environmental quality and for 
sustainable economic development in long run. The sectoral 
analysis of EKC would be helpful for designing a comprehensive 
growth, energy and environmental policy to maintain living 
standard of Turkish people. 
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